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As a step in the direction of "computer vision," several programs 
have been written which transform the output of a vidissector into some 
mathematical descriptions of the boundaries enclosing the objects in the 
field of view. Most of the discussion concerns the techniques used to 
transform a sequence of points, presumably representing a curve in the 
two-dimensional plane of view, into the best-fit conic-curve segment, or 
best-fit straight line. The resultant output of this stage is a list of 
such segments, one list for each boundary found. 
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As a step in the direction of "computer vision/’ several programs 
have been written which, applied successively, transform the point-view 
from the vidissector into a mathematical description of the boundaries 
of the objects in the field of view. An executive program will co-ordinate 
the actions of the various programs described below in order to find each 
and every object in the view. For example, several attempts may have to 
be made in order to be sure that a single object has been seen, using such 
clues as color contiguity, graininess, or texture, light intensity, etc. 

The programs REGI0NS1 and BNDS0RT1 examine the points in the current 
subview of the vidissector and determine which ones belong to the object 
or region being sought. The executive determines which section of the 
vidissector to look at, the density of points to be interviewed, and the 
predicates that define a region. The end result is an ordered list of 

points consituting the boundary of the region. (Each point is a list of 
two LISP numbers, which are the x- and y - coordinates on the vidissector 
screen of the geometric point represented.) Further information on these 
routines may be obtained from Gerald Sussman. 

POLYSEG takes the output of BNDS0RT1 and a straight-line approxima¬ 
tion to the boundary of the figure, using techniques adapted from the work 
of Roberts (see Lincoln Laboratory Technical Report No. 315, "Machine Per¬ 
ception of Three-Dimensional Solids"). The purpose of this routine is 
twofold: (1) to smooth out minor pieces of noise in the data, and (2) to 

reduce the total number of points representing the boundary. Hopefully 
the error tolerance will insure that a sufficient number of short straight 
lines are used to approximate a curved line segment on the boundary. 
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Perhaps more work could be put into this routine in the form of a better 
heuristic determination of corners as they appear in a boundary. POLYSEG 
returns a list of n+1 points for an n-sided polygonal approximation to 
the boundary of the region. The first point of the list is duplicated 
as the last point also, thus explaining the additional point. 

ANGLES takes the output of POLYSEG and returns a list of four items 

(1) Augmented polygon list [see below] , 

(2) Reasonable approximation to the diameter of the region described 
by the input polygon list, 

(3) "CLCK" or "CNTK" indicating clockwise or counterclockwise ordering 
of the boundary points, 

(4) Area of the approximating polygon (positive,.floating point). 

The original polygon list, i.e. the input list, has only the x-and y-coor- 
dinates for each vertex; the augmented list,(1) above, provides also the 
inside angle at each vertex and the distance to the next-in-order vertex. 
Thus each member (x^y_^) °f the input becomes (x^y_^a_^d^) t ^ ie augmented 
output list. 

SLPCS takes the output of ANGLES (providing the points are ordered 
clockwise) and attempts to segment the boundary into true straight lines 
or curved pieces. By "piece" we mean a segment of the boundary which, 
according to the following criteria, is a straight line, or a straight- 
line approximation to some curve without inflections. The output is then 
a listing of the pieces of the boundary in order, and the listing of the 
points of each peice is ordered the same as in the input. In addition to 



the points of the relevant boundary segment, the first member of each 
piece (which is otherwise a list of points) is a fixed point integer 0, 
1, or 2 corresponding to "straight line," "curved out," and "curved in" 


respectively. Example: 




Notice that each piece is complete in itself so that its endpoints are 
the same points as will be found on one end of either adjoining piece. 
Several heuristics are applied to decide where one piece should end and 
another begin. Clearly an inflection point will be a breakpoint. Some 
thresholds are computed: LNTSH1 is such that at any point on the boundary 
if the distance to the next point is greater than LNTSH1, then this seg¬ 
ment is automatically taken to be a straight line piece. If the said dis- 

tnace is less than LNTSH2 then it is assumed that this will not be a 

straight line piece, but rather part of some curve piece. The other cases 

are resolved by some additional formulae: consider point i with associ¬ 
ated distance d^ and inside angle (outside angle a_^ = tt - a^) . 
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We will consider the segment (x.y.) 

piece if d. > LNTSH1 or if 

x 


(x i+i y i+i ) 


to be a straight line 


[d. > LNTSH2 /\ (a. n ' d. > CRVF v a? • d. > CRVB)]. 

1 1+1 1 11 

Any curved piece is extended until an inflection point or straight line 
is reached. From the above formula, it is clear that straight lines are 
usually found by the local property of "sharp corner”. The formula involves 
the two endpoints of a straight line in slightly different ways because it 
was found to be efficacious in discerning certain odd types of corners while 
using a one-way scan of the points from ANGLES. For example 



Experiments were carried out to determine suitable values for the above 
thresholds, with the following results: 
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LNTSH1 = 0.25 * DIAMETER 
LNTSH2 = 0.03 * DIAMETER 
CRVB = 0.022 * DIAMETER 
CRVF = 0.02068 * DIAMETER 

The curved pieces output by SLPCS may be represented as some segment 
of a conic (general second degree curve) since there are no inflection 
points in such curves (some pathological cases do arise however, such as 
egg—shaped closed curves and spiral sections, but for the moment these 
are ignored, and probably will not be parts of viewed objects anyway). 
BSTCRV accepts any piece from the output of SLPCS, including straight 
lines which it leaves unaltered, and returns a list of eight members repre 
seating a conic-curve segment. That is, for straight lines 

(0 (x lYl ) (x 2 y 2 )) BSTCRV^ (0 (xiyi) (x 2 y 2 )) 

and for curves which can be approximated by a circular or,elliptical arc 
(IND (x!y x ) ... (x 2 y 2 )) BSTCRV. (IND 'ELLPSE' P 0 A P'iS) 

where 


IND 

is 

1 or 2 for curved out or in 


’ELLPSE' 

is 

the type of curve, as opposed to 

hyberbolic 

Po 

is 

a point representing the center 

of the ellipse 

1 — ' 

is 

a LISP number for the length of 

the principal 



semi-axis 


/2 

is 

for the secondary semi-axis 
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ot is an angle -tt/2 < a tt/2 measuring the tilt of the 

principal axis from the horizontal, about the center 
of the ellipse. 

Pi, 3 specify the segment of the curve in the following 

manner: "beginning at point Pi and using Pq as cen¬ 

ter point, take B radians about Pq, in the clock¬ 
wise direction if IND = 1 otherwise in the counter¬ 
clockwise direction. This is the relevant segment. 


Clearly for a circle x_i = X 29 ami for any complete circle or ellipse, 
3 = 2tt and Pi is any point on the curve. 

Example 



-„- > 

x 

For second-degree curves which are not ellipses, i.e. hyperbolae or 
parabolae we simply return the general-equation coefficients and the 
first and last point on the arc. 


(IND (x,vi) ... (x v 1) BSTCRVv (IND 'HYPER' (A B) C D E (x 1 v 1 ') (x v ')') 

n n * v x 1J1 nn 

a 

where the best fit curve is 

Ax 2 + Bxy + Cy 2 +Dx + Ey = 1 
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Appendix I 

Atoms Used as Free Variables 
Throughout the Vision Programs 


PI ) Set to it and 2 tt by a top level call to SETQ, these 

TWOPI ) atoms provide convenient reference to some familiar 

numerical values. 


CTOL 

ETOL 


BSTCRV requires the mean'square error in an approximation 
to be less than some tolerance. Since circles are intrin¬ 
sically more simple than general conics, a circular fit is 
first examined: if the computed circle has center (a,b) 
and radius v, then it is accepted iff 






2 

- (1 + 


a 2 + b^ 
“72 - 



< CTOL 


If the circle is not accepted, a general second-degree 
curve fit is tried and accepted iff 


1 

n 


n 


i=l 


Ax t + Bx.y. 4- Cy. + Dx. + Ey. -1 

l li l l ^ i 


< ETOL 


See Appendix II for the formulae by which the best fit 
curve coefficients are determined. 

Suitable values for these tolerances are ETOL = .01 

CTOL * .0015 

but some applications may want to vary them. 



} 


i 


TOL 


HI 


N2 


POLYS1G determines a straight lias approximation to a 
curva by initially *l«|$ag up Klpeints and extending 
by N2 points at a tine until rhu ipnnel tied sean-s^uare 
error of the beat-fit strAAght line to these points is 

greater than T*L. "lggT~yiT w meant least-squareg■best 

• ■ 

fit in which 


i ® 

is 

»£1 


(ax. + by 



appropriate choice 


of a,b, and c. 
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Appendix II 

BEST FIT CIRCLE 

For a data set ^Cx^y^) ... (x y )J the best fit circle, with equa¬ 
tion ax 2 - 4- bx 4- ay 2 + cy = 1.0, is the one which minimizes 

2 

ax? + bx. + ay? + cy. - 1.0 
l l l 

Differentiating the error expression above with respect to a, b, and c, 
leads to three simultaneous linear equations, which then can easily be 
solved. The center of the circle is at 

•j^ j -T - 

(- - tt") and it has radius \ / — + A 2 + B 2 

2a 2a V ,a 



BEST FIT CONIC 

As above, we solve five simultanious linear equations in order to 
minimize 



+ bxy + cy 


2 


+ dx + ey -1.0] 2 . 


Actually, a parametric form of the second-degree equation would be more 
revealing, 

ax 2 + bxy + cy 2 + dx + ey = f. 

One must place some additional constraint, however, since there*are only 
five degrees of freedom in a conic curve, but the parametric equation has 
six unknowns. Imposing 

\l a 2 + b 2 + c 2 + d 2 + e 2 + f 2 = 1 


looks plausible, but the resulting simultaneous equations are horribly 
non-linear, so the much simpler condition f = 1 has been chosen. Cases in 
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which this latter assumption leads to numerical nonsense rarely occur 
in practice. 


DETERMINATION OF ELLIPSE PARAMETERS 

If b - 4ac < 0, then the curve represented by 

•) o 

ax 4- bxy + cy + dx + ey = 1.0 

is an ellipse. Otherwise it is a parabola or hyperbola. If the whole 
plane is rotated about the origin through an angle of -a radians, the 
axes of the ellipse will be parallel to the x- and y-axes. Thus consider 
the curve 

Ax 2 + Cy 2 + Dx + Ey = 1.0 


where 


a = sign (b) 


sign(a-c) 


1 

2 tan 


-1 



a-c 


sign(x) = 


1 if x > 0 
I 

-1 if x < 0 


o 9 

A = a cos ot 4- b cosot sinot + c sin ot 

9 9 

B = b(cos ot - sin ot) + 2(c-a) sinot cosot = 0 

2 9 

C = a*sin ot - b*cosot • sinot 4- c - cos^ot 
D = d cosot + e sinot 
E = -d sinot + e cosot 

This represents an ellipse with center at (- 

M 1JL 

semi-axis length is 



whose horizontal 




+ d 2 c 2 


+ ACE 2 


2AC 


and whose vertical semi-axis 


length is 
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